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Abstract 

Recent observations of Galactic Black Hole Candidates (BHCs) suggest that those that 
are superluminal jet sources have more rapid black hole spin rates than otherwise normal 
BHCs. This provides observational support for models of relativistic jet formation that extract 
rotational energy of the central black hole. To investigate this mechanism, we have developed 
a new general relativistic magnetohydrodynamic code in Kerr geometry. Here we report on the 
first numerical simulation of the formation of a relativistic jet in a rapidly-rotating (a = 0.95) 
Kerr black hole magnetosphere. We assume that the initial velocity of the disk is zero. We 
find that the maximum velocity of the jet reaches 0.93c (Lorentz factor, 2.7) and the terminal 
velocity of the jet is 0.85c (Lorentz factor, 1.9). On the other hand, for a non-rotating (a = 0) 
Schwarzschild black hole, the maximum jet velocity is less than 0.6c as far as we calculated with 
the similar condition of the magnetosphere to that of the Kerr black hole case. These numerical 
results show the importance of the rapidly rotating black hole for the relativistic jet formation 
clearly. 



1 Introduction 

Radio jets ejected from radio loud active galactic nuclei (AGNs) sometimes show proper motion 
with apparent velocity exceeding the speed of light c fil Ej. The widely-accepted explanation 
for this phenomenon, called superluminal motion, is relativistic jet flow in a direction nearly 
along the observer's line-of-sight with a Lorentz factor greater than 2.0 Q. Such relativistic 
motion is thought to originate from a region very close to the putative supermassive black hole 
which is thought to power each AGN g]. On the other hand, the great majority of AGN are 
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radio quiet and do not produce powerful relativistic radio jets 0|. These two classes of active 
objects (radio loud and quiet) are also found in the black hole candidates (BHCs) in our own 
Galaxy. Objects with superluminal jets, such as GRS 1915+105 and GRO J1655-40, belong to 
the radio loud class ](| [?|. Other objects such as Cyg X-l and GS 1124-68 are relatively radio 
quiet and produce little or no jet. 

What causes the difference between the two classes? Recent observations of the BHCs in 
our Galaxy suggest that the Galactic superluminal sources contain very rapidly rotating black 
holes (specific angular momentum of rotating black hole a = 0.9 — 0.95), while the black holes 
in Cyg X-l and GS 1124-68 are spinning much less rapidly (a — 0.3 — 0.5) ||. According to 
recent (nonrelativistic) studies of magnetically-driven jets from the accretion disk by Kudoh & 
Shibata j^, [u| , the terminal velocity of the formed jet is comparable to the rotational velocity of 
the disk at the foot of the jet. Further nonrelativistic simulations of jet formation confirm these 
results |ll], [l^], except for the extremely large magnetic field case fl5| j in which very fast jets 
can be produced. The rotation velocity at the innermost stable orbit of the Schwarzschild black 
hole (r = 3rs) is 0.5c. (Here rs = 2GA1bh/c 2 is the Schwarzschild radius, where G and Mbh 
are the gravitational constant and black hole mass, respectively.) In addition, recent studies 
indicate that the poloidal magnetic field strengths in disks around non-rotating black holes may 
be rather weak [tt4j]. Therefore, a jet produced by MHD acceleration from an accretion disk 
around a non-rotating black hole should be sub-relativistic and very weak. In fact, numerical 
simulations of jet formation in a Schwarzschild metric show only sub-relativistic jet flow Jl5| |, 
except for the case with the artificial condition of a hydrostatic corona in equilibrium |r| ]. 

Several mechanisms for relativistic jet formation from rotating black holes have been pro- 
posed ^o) . However, up until now no one has performed a self-consistent numerical 
simulation of the dynamic process of jet formation in a rotating black hole magnetosphere. 
To this end, we have developed a Kerr general relativistic magnetohydrodynamic (KGRMHD) 
code. In this paper we report briefly on what we believe is the first calculation of its kind 
— simulation of relativistic jet formation in a rotating black hole magnetosphere. When the 
rotation parameter of the rotating black hole a = 0.95, the maximum velocity of the formed 
jet reaches 0.93c (Lorentz factor, 7 = 2.7), and the terminal velocity is 0.85c (7 = 1.9). On the 
other hand, when a vanishes, only a slow jet is formed even if the initial condition of magne- 
tosphere is almost the same as that of the Kerr black hole case. This result clearly shows the 
importance of the rapidly spinning black hole for the formation of relativistic jet. 



2 Numerical Method 

We use a 3 + 1 formalism of the general relativistic conservation laws of particle number, 
momentum, and energy and Maxwell equations with infinite electric conductivity [^lj . The Kerr 
metric, which provides the spacetime around a rotating black hole, is used in the calculation. 
When we use Boyer-Lindquist coordinates, x° — ct, x 1 = r, x 2 — 6, and x 3 = 0, the Kerr 
metric g M „ is written as follows, 
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ds 2 = g^dx^ dx v — — /ig(cdt) 2 + S ' h 2 (dx 1 ) 2 — 2h^2,cdtdx 3 . (1) 



By modifying the lapse function in our Schwarzschild black hole code (a = ^Jl — rs Jr) to 
be a = y/h,Q + f^, and adding some terms to Q3, we were able to develop a KGRMHD code 
relatively easily. (See Appendix C in jlB) for more details on this procedure and the meaning 
of symbols used.) 

We use the Zero Angular Momentum Observer (ZAMO) system for the 3- vector quantities, 
such as velocity v, magnetic field B, and so on. For scalars, we use the frame comoving with 
the fluid flow. The simulation of the Kerr (Schwarzschild) black hole is performed in the region 
0.75rs (1.2rs) < r < 20rs, < 9 < tt/2 with 210 x 70 mesh points, assuming axisymmetry 
with respect to the z-axis and mirror symmetry with respect to the plane z — 0. A simplified 
radiative boundary condition is employed at r = 0.75rs (r = 1.2rs) and r — 20rs for the 
Kerr (Schwarzschild) black hole case. In the simulations, we use simplified tortoise coordinates, 
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Figure 1: A relativistic jet formation around a Kerr black hole (rotation parameter a — 0.95). (a) 
Initial, (b) Final stage at t = 65rg. Color shows the logarithm of mass density. Vector indicates 
the velocity. The lines show the magnetic field. The broken line near the horizon shows the inner 
boundary of the calculation region. The relativistic jet are formed at t = 65rs. The maximum 
velocity of the jet is 0.93c (Lorentz factor, 7 = 2.7) and the terminal velocity is 0.85c (7 = 1.9). 
The solid straight line shows the plot line used in Fig. 3. 

x = log(r/rH — 1), where ru is the radius of the black hole horizon. Here rrr = 0.656rs 
(m = rs) when a — 0.95 (a = 0). To avoid numerical oscillations, we use a simplified TVD 
method Jl^, Q . We checked the KGRMHD code by computing Keplerian motion around 

a rotating black hole and comparing with analytic results. 

3 Results 

Figure 1 shows the initial and final states of the calculation of relativistic jet formation near 
a Kerr black hole. These figures show the rest mass density (color), velocity (vector), and 
magnetic field (solid lines) in < R < 8rs, < z < 8rs- The poloidal magnetic field lines 
are drawn as the contour lines of the azimuthal (<j>) component of the four vector potential 
(A t , A r , Ag, A$). The black region at the origin shows the black hole horizon. The specific 
angular momentum of the black hole is a = 0.95 and the radius is ru = 0.656rs. The initial 
state in the simulation consists of a hot corona and a cold accretion disk around the black hole 
(Fig. la). In the corona, plasma are assumed to be nearly stationary falling with the specific 
enthalpy h/pc 2 = 1 + Tp/[(T — l)pc 2 ] = 1.3, where p is the rest mass density, p is the pressure, 
and r is specific heat ratio and set T = 5/3. Far from the hole, it becomes the stationary 
transonic solution exactly. The accretion disk is located at |cot#| < 0.125, r > ro = 3rs and 
the initial velocity of the disk is assumed to be zero. The mass density of the disk is 300 times 
that of the corona. In addition, the magnetic field crosses the accretion disk perpendicularly. 
We use the vector potential A^ of the Wald solution to set the magnetic field, which 
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Figure 2: A sub-relativistic jet formation around Schwarzschild black hole (a = 0). (a) Initial (b) 
Final stage at t = 65rs. Color shows the logarithm of mass density. Vector indicates the velocity. 
The lines show the magnetic field. The broken line near the horizon shows the inner boundary of 
the calculation region. The maximum and terminal velocities of the jet are 0.6c and 0.45c, both 
sub-relativistic, respectively. The solid straight line shows the plot line used in Fig. 4. 

produces a uniform magnetic field far from the Kerr black hole. However, we do not use At 
from the Wald solution; instead, we use the ideal MHD condition E + vxB = 0to determine 
the electric field E. Here the Alfven velocity and plasma beta value at the disk (r = 3rs) are 
wa = 0.01c and /3 ~ 12, respectively. 

Figure lb shows the state at t = 65rs, where rs is defined as rs = rs/c. The accretion disk 
falls into the black hole rapidly, and the disk plasma crosses the horizon, as suggested by the 
crowded magnetic field lines near the horizon. The magnetic field lines become radial due to 
the dragging by the disk near the black hole. The jet is formed almost along the magnetic field 
lines, while outside of the jet, the magnetic field lines are bent strongly. The maximum velocity 
of the jet reaches to 0.93c (Lorentz factor, 7 = 2.7) at R = 5.3rs, z — 4.4rs and its terminal 
velocity is 0.85c (7 = 1.9). The terminal velocity is defined by the jet velocity far from the 
black hole (r ~ 10rs). We found the jet has a two-layered-shell structure with respect to the 
density. The low density part of the jet is fast and relativistic, while the high density part at 
the outer shell is slower and sub-relativistic. The fast jet dominates the jet and the slow, dense 
jet is thin. The high density jet is connected to the high density region of the accretion disk. 

To show the importance of the black hole rotation for relativistic jet formation, we also 
performed a simulation of the Schwarzschild black hole case (a = 0) (see Fig. 2). The initial 
condition is almost the same as that of Kerr black hole case except for the black hole rotation 
(Fig. 2a). In the initial stage, the coronal plasma is assumed to be stationary falling with 
h/pc 2 — 1.3 (see Fig. 2a). The disk is located at the same place of the Kerr black hole case: 
|cot#| < 0.125, r > ro = 3rs, and the density of the disk is also 300 times larger than that of 
corona. The uniform magnetic field crosses the disk perpendicularly. The strength is the same 
as that of the Kerr case. Figure 2b shows the state at t = 65rs . The maximum velocity of the 
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jet is 0.6c and the jet terminal velocity is 0.45c, both sub-relativistic. The jet also has two- 
layered-shell structure. The inner shell has low density, while the outer shell has high density. 
The high density jet dominates the jet, which is different from the Kerr black hole case. The 
disk plasma also falls and crosses the event horizon rapidly as the Kerr black hole case. Note 
that the inner edge of the high density region is not identical with that of disk plasma. The 
plasma of the disk passes through the edge and falls into the black hole. The small high pressure 
region is formed at R ~ 3rs in the disk. The high pressure is caused by the shock in front of 
the region. The dense jet is accelerated by the high pressure. This mechanism is the same as 
that of the case of the Kepler disk around the Schwarzschild black hole |l6j . The magnetic field 
lines are also bent by the jet strongly. 

To show the structure of the relativistic jet in the Kerr black hole case, various physical 
quantities on the line z — R — OArs (1.3rs < R < 6.3rs) are shown in Fig. 3. The line is almost 
along the fast, low density jet (see straight line in Fig. lb). Figure 3a shows the magnetic 
pressure from the toroidal component of the magnetic field, £>^/2p is very large near the black 
hole (1.3rs < R < 2.5rs) where p; is the magnetic permeability of the vacuum. Because initial 
magnetic field is = 0, B z = 0.15poc\/JI, where po is the unit mass density, the field increases 
to more than 60 times the initial magnetic field. This amplification is caused by the shear of 
the plasma flow observed by the Boyer-Lindquist frame or local nonrotating frame (LNRF) due 
to the frame dragging effect of the rota ting black hole . This effect is called a metric- 

shear-driven MHD dynamo by Meier p8|. With the assumption that the azimuthal velocity 
component is zero, general relativistic Faraday law of induction and ideal MHD condition yield, 

^f=faB r + faB e , (2) 

where /i = c(h3/hi)d(Q,3/h3)/dr, and fa = c(/i3 / '}i2)d(Q,3 / '/13) / '88 . Note that fa is one order 
smaller than fa when a ~ 1. The ridge of the fa profile is located near the spherical surface, 
r — rs and the magnetic pressure has a maximum near the surface. Where does the energy of the 
amplification of the magnetic field comes from? The energy does not come from the gravitational 
energy or thermal energy of the disk, because this effect occurs even when the plasmas of the 
disk and corona are rest and cool. The only other possible energy source is the rotation of the 
black hole itself. In fact, the increase in the azimuthal magnetic field component depends on 
the shear of the rotational variable SI3. We conclude that the energy of the amplification of 
the magnetic field is supplied by the extraction of the rotational energy of the black hole. Note 
that this mechanism of the extraction of the rotational energy is different from the previous 
electromagnetic mechanisms jl^, |l9| . The detail comparison with these models will be presented 
in our following paper. If there is no disk, the magnetic tension causes the azimuthal motion 
of the plasma to prevent the dynamo. The heavy disk plays an important role to keep the 
magnetic energy from the energy of the black hole rotation to produce a relativistic jet. We 
estimate the power from the electro-magnetic field, Wem = v-(E + JxB) and the gas pressure, 
W sp — — v- Vp (Fig. 3c). The power contributions of the electromagnetic field and gas pressure 
are comparable. The high gas pressure is caused by the shock in the disk at r ~ 2rs. The 
density and the pressure are small in the fast jet region (3.5rs < R), which means that the 
jet comes from the corona, not the disk. Both the strong magnetic and gas pressures push the 
low density plasma to become the relativistic jet. Figure 3b shows the azimuthal velocity is 
small over the whole range. This indicates the jet acceleration is caused by the magnetic and 
gas pressure rather than the magnetic tension (or centrifugal force). The magnetically driven 
part of the mechanism is the same as that of Uchida & Shibta p9j and Shibata & Uchida 
]. Punsly & Coroniti (1990) proposed an alternative magnetic force driven mechanism of the 
relativistic jet formation. They considered the magnetic tension causes the azimuthal velocity 
to produce the jet due to centrifugal force like the model of non-relativistic jet formation from 
the accretion disk |3lJ. Figure 3b shows that in the fast jet region (4rs < R), the azimuthal 
velocity is small but positive (v^, > 0), which is remnant of the frame dragging by the rotating 
black hole. In the acceleration region (R < 4rs), on the other hand, the azimuthal velocity is 
negative, as predicted by Meier ||2fi}| . Accretion of this negative angular momentum material 
ultimately reduces the spin of the black hole and is at least partially responsible for extracting 
the hole's rotational energy. 
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Figure 3: Physical quantities of the Kerr black hole case along z = R— 0.4rs (1.3rs < R < 6.3rs) 
at t = 65rs (see straight line in Fig. lb), (a) Mass density, p, gas pressure, p, and poloidal magnetic 
pressure, S|/2. (b) Poloidal velocity, v p and azimuthal velocity v<f,. (c) Power of electromagnetic 
field (Wem) and gas pressure (W gp ) to the acceleration of the jet. 
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Figure 4: Physical quantities of the Schwarzschild black hole case along z = R — 2rs (2.3rs < 
R < 7.3rs) at t = 65rs (see straight line in Fig. 2b). (a) Mass density, p, gas pressure, p, and 
poloidal magnetic pressure, Bi/2. (b) Poloidal velocity, v p and azimuthal velocity, v^. (c) Power of 
electromagnetic field (Wem) and gas pressure (W gp ) to the acceleration of the jet. 
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Figure 4 shows the physical variables of the Schwarzschild case at t = 65rs on the line 
z — R — 2rs (2.3rg < R < 7.3rg). The line is also almost along the main, dense jet (see straight 
line in Fig. 2b). The terminal velocity of the jet is 0.45c, sub-relativistic (Fig. 4b). However, 
the mass density and pressure are relatively high (Fig. 4a). This means that the jet comes 
from the disk. The jet is accelerated only by gas pressure (Fig. 4c). The high pressure is 
caused by the shock heating in the disk as the same as in the case of a Keplerian disk around 
a Schwarzschild black hole llq] . It is reasonable that the azimuthal components of velocity and 
magnetic field vanish because both the black hole and disk have no rotation. In such a case, 
magnetic acceleration is not permitted and only a gas pressure driven jet may be formed. 

4 Discussion 

We have also performed a maximally rotating black hole case (a = 1). Preliminary results show 
that a strong relativistic jet is formed, with a maximum velocity of 0.98c (Lorentz factor 5!) 
and a terminal velocity of 0.9c (Lorentz factor, 2.3). Unfortunately, the inner boundary of this 
calculation is located at r — 0.8rs, which is a little bit far from the horizon, th = 0.5rs. We 
believe that a faster relativistic jet will be formed when we use an inner boundary nearer to 
the horizon. On the other hand, in the Schwarzschild case with similar initial conditions in 
the magnetosphere, only a sub-relativistic jet is formed. We have also performed a number of 
simulations with rotating disks around Schwarzschild black holes jl5| . In these cases, all formed 
jets are sub-relativistic, except for a special case with hydrostatic equilibrium corona [jl6[] . 

Of course, more simulations of jets formed by rotating black holes accreting magnetized 
plasma need to be done in order to fully understand this formation mechanism. However, our 
numerical results show clearly that the spin of the black hole is important in the formation a 
relativistic jet. A comprehensive model of relativistic jet formation around Kerr black hole, 
along with more detailed data analysis, will be presented in our following paper. 

Recent X-ray observations of black hole candidates (BHCs) in our Galaxy show the relation- 
ship between the jet speed and the rotation of the central objects || . A quasi-periodic-oscillation 
(QPO) is found in the power spectrum of the Galactic superluminal source GRO J1655-40 [ [j2"[ 
and GRS 1915+105 133] . Using QPO observations, it is concluded that the both Galactic super- 
luminal sources contain very rapidly rotating black holes (a ~ 0.95). Similarly, using the QPOs 
in Cyg X-l 0] and in GS 1124-68 (35), it is shown that the black holes in these two sources 
are spinning less rapidly (a ~ 0.35 and a ~ 0.48, respectively). Based on these results, Cui, 
Zhang, & Chen (1998) proposed that the difference between BHCs that are also superluminal 
jet sources and otherwise normal BHCs is in the spin of their black holes. This observational 
result is consistent with our numerical results. 
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manuscript. One of author (S. K.) thanks M. Inda-Koide for her discussion and important 
comments for this study. We also thank M. Takahashi, A. Tomimatsu, J. F. Hawley, R. D. 
Blandford, M. C. Begelman, P. E. Hardee, J.-I. Sakai, and R. L. Mutel for their discussions and 
encouragement. We appreciate the support of the National Institute for Fusion Science and the 
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